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Abstract. The present study aimed to investigate the gene 
functions and expression profiles in perihematomal (PH) 
brain regions following spontaneous intracerebral hemor- 
rhage. The gene expression profiles were downloaded from 
the Gene Expression Omnibus database under accession 
number GSE24265, which includes 11 brain samples from 
different regions, including four samples from PH areas, four 
from contralateral grey matter (CG) and three from contra- 
lateral white matter (CW). The gene expression profiles were 
pre-processed and the differentially expressed genes (DEGs) 
between PH and CG tissue, and PH and CW tissue were identi- 
fied using R packages. The expression of genes in different 
tissues was analyzed by hierarchical clustering. Then, the 
interaction network between the DEGs was constructed using 
String software. Finally, Gene Ontology was performed and 
pathway analysis was conducted using FuncAssociate and 
Expression Analysis Systematic Explorer to identify the gene 
function. As a result, 399 DEGs were obtained between PH 
and CG, and 756 DEGs were identified between PH and CW. 
There were 35 common DEGs between the two groups. These 
DEGs may be involved in PH edema by regulating the calcium 
signaling pathway [calcium channel, voltage-dependent, 
T-type, all subunit, Ca^Vcalmodulin-dependent protein 
kinase II a (CAMK2A), ryanodine receptor 2 (RYR2) and 
inositol 1,4,5-trisphosphate receptor, type 1 (ITPRl)], cell 
proliferation (sphingosine kinase 1), neuron differentiation 
(Ephrin-A5) or extracellular matrix-receptor interaction 
[collagen, type I, a 2, laminin Bl (LAMBl), syndecan 2, 
fibronectin 1 and integrin a5 (ITGA5)] . A number of genes 
may cooperate to participate in the same pathway, such as 
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ITPR1-RYR2, CAMK2A-RYR2 and ITGA5-LAMB1 interac- 
tion pairs. The present study provides several potential targets 
to decrease hematoma expansion and alleviate neuronal cell 
death following spontaneous intracerebral hemorrhage. 

Introduction 

Spontaneous intracerebral hemorrhage (ICH) involves the 
extravasation of blood within the brain parenchyma in the 
absence of trauma or surgery, accounting for 10-30% of all 
strokes worldwide (1). ICH is a complex, dynamic process 
that consists of three distinct phases (2): i) Initial hemorrhage, 
retraction of the clot begins in the first few hours following 
ICH. ii) Hematoma expansion, hematoma expansion occurs 
24 or 48 h following the initial onset of symptoms, which 
involves a sudden increase in intracranial pressure, disrupting 
the integrity of the local tissue and blood-brain barrier. In 
addition, obstructed venous outflow promotes the activation 
of tissue thromboplastin, resulting in local coagulopathy, 
(iii) Perihematoma (PH) edema, PH cerebral edema develops 
within the first few days following ICH as a result of an 
inflammatory response [activation of cellular (leukocytes, 
macrophages and microglia) and molecular (cytokines, prote- 
ases and reactive oxygen species) components] (3) secondary 
to local release of thrombin and other end products of coagula- 
tion from the hematoma, and also due to cytotoxic mediators 
(hemoglobin and its degradation products, heme and iron) that 
are released from red blood cells (4-6). 

The hematoma expansion rate and PH edema volume are 
key factors to predict neurological deterioration and survival. 
Hematoma volume >30 ml is associated with increased 
mortality (7). When the hemorrhagic volume exceeds 150 ml, 
cerebral perfusion pressure falls to zero, which subsequently 
leads to patient mortality (8). Therefore, prevention of hema- 
toma expansion is a logical therapeutic target (9). Additionally, 
it is possible that anti-inflammatory and anti-apoptotic drugs 
may have the therapeutic potential to ameliorate secondary 
brain injury following ICH. For example, sesamin confers 
neuroprotection in rat ICH by suppressing activation of 
microglia and the p44/42 MARK pathway (10). Ginkgolide B 
was able to significantly suppress the gene expression of 
Toll-like receptor 4 and nuclear factor-KB, reduce the concen- 
trations of tumor necrosis factor-a, interleukin (IL)-ip and 
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IL-6 as well as reduce the number of apoptotic neuronal cells 
in hemorrhagic rat brain tissues (11). 

Recently, Carmichael et al (12) attempted to investigate the 
gene expression pattern in PH tissues of ICH rat models in a 
microarray study, which simultaneously measured transcript 
abundances of thousands of genes in a cell population or 
tissue. The results indicated a robust upregulation of pro-/anti- 
inflammatory networks and a downregulation of neuronal 
signaling pathways (12). Furthermore, Rosell et al (13) investi- 
gated the genomic profile following spontaneous human ICH 
and suggested that the significantly overexpressed genes in 
the PH areas encoded for cytokines, chemokines, coagulation, 
cell growth and proliferation factors, while the downregulated 
genes encoded proteins involved in the cell cycle or neuro- 
trophins. Therefore, the identified genes may be potential 
targets for ICH therapy. 

The present study aimed to further identify the differen- 
tially expressed genes (DEGs) between PH and contralateral 
healthy tissues [grey (CG) matter and white (CW) matter] 
using a different statistical approach and threshold. Further 
interaction network and function analyses were also performed 
for these DEGs. It is expected that these data may further 
elucidate the processes and mechanisms of death caused by 
hematoma following spontaneous ICH. 

Materials and methods 

Microarray data. The gene expression profiles were down- 
loaded from the Gene Expression Omnibus under accession 
number GSE24265 (13), which includes 11 brain samples from 
different regions of four deceased patients who had suffered 
from a supratentorial ICH within four days prior to mortality. 
Of these samples, four were from PH tissues, four from CG 
and three from CW. The study group included two female 
and two male subjects with a median age of 79 years (range, 
68-92 years). The samples from PH, CG and CW tissue were 
obtained within the first 5 h following mortality, snap frozen 
in liquid nitrogen and stored at -80°C until RNA isolation. The 
original studies by Rosell et al (13) was approved by the Ethics 
Committee of the Vail d'Hebron Hospital (Barcelona, Spain) 
and informed and written consent was acquired from all of 
the patients or patient's relatives. Demographic and tissue 
sampling data are presented in Table I. The array platform was 
GPL570 [Human Genome U133 Plus 2.0 Array (Affymetrix 
Inc., Santa Clara, CA, USA)]. The array annotation informa- 
tion was from Affymatrix, including the information for all 
probes on Affymetrix ATHl (25K; (Affymetrix Inc.). 

Data pre-processing and identifying DEGs. Firstly, the raw 
data were transformed into a readable expression profile 
format and then the missing data were imputed according to 
Troyanskaya's method (14) followed by normalization (15). 
DEGs between PH and CG, and PH and CW were identified 
using the limma package in R (16). Benjamini-Hochberg 
(BH) multiple test correction was then used to adjust the 
P-values (17). Only the genes with a false discovery rate (FDR) 
<0.05 and llogFC(fold change)! >1 were selected as DEGs. 

Tissue specific DEGs. The expression of genes is time- and 
space-specific, and in different tissues may be significantly 



Table I. Demographic and tissue sampling data. 



Case 


Sex 


Age 


Time S-E 


Time E-T 


PH 


CG 


CW 


27 


F 


72 


12h 


5 h 


1 


1 


1 


28 


F 


95 


39 h 


2.5 h 


1 


1 


1 


30 


M 


67 


17 h 


3 h 


1 


1 




31 


M 


86 


4 days 


4h 


1 


1 


1 



F, female; M, male; S-E, stroke to exitus; E-T, exitus to tissue samples; 
PH, perlhematomal ai'eas; CG, contralateral grey matter; CW, contralat- 
eral white matter. 



different. The genes were clustered based on their expression 
in two groups using Cluster analysis (18). The expression 
differences in the various tissues were visualized by heatmap 
(Clifton Watt, Stanford University, USA). 

Construction of the interaction network. Numerous activi- 
ties in the body, including important physiological activities 
and responses to the external or internal environment, are 
performed by the binding and dissociation of proteins. All 
of these activities are conducted by signaling transduction 
networks consisting of protein-protein interactions (19). 
Therefore, detailed investigations into protein interac- 
tion networks is essential for the understanding of bodily 
systems (20). Based on the sequence and structure informa- 
tion of DEGs, their interaction possibility was predicted 
according to the String database (http;//string.embl.de/) that 
quantitatively integrates interaction data from four sources 
[genomic context, high-throughput experiments, (conserved) 
co-expression and previous knowledge] for a large number of 
organisms, and transfers information between these organ- 
isms where applicable (21). Furthermore, the interaction 
network was also constructed for the two groups; PH vs. CG 
and PH vs. CW. 

Functional enrichment of DEGs in interaction networks. 
The annotation was performed for the genes in protein inter- 
action networks. FuncAssociate (http://llama.med. harvard, 
edu/funcassociate) is a web application that identifies proper- 
ties enriched in lists of genes or proteins that emerge from 
large-scale experimentation (22). The interactive DEGs were 
screened and functional enrichment analysis was performed 
using FuncAssociate (22). FDR<0.05 was set as the threshold 
for this analysis using the hypergeometric distribution. 

Pathway analysis of DEGs in interaction networks . Expression 
Analysis Systematic Explorer (EASE, http://david.niaid.nih. 
gov/david/ease.htm) is a customizable, standalone software 
application that facilitates the biological interpretation of 
gene lists derived from the results of microarray, proteomic 
and Serial Analysis of Gene Expression studies (23). EASE 
queries the Kyoto Encyclopedia of Genes and Genomes 
pathways database (http://www.genome.ad.jp/kegg/) and 
determines functional enrichment by calculating a Fisher's 
exact test P-value for each pathway. P<0.05 was considered to 
indicate as significant gene- enrichment in a specific annota- 
tion category. 
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Figure 1. Box plots of normalized gene expression data. Blue, pink, orange, 
respectively represents brain PH, CG and CW. The median of normalized 
data are almost at the same level, indicating good standardization. PH, peri- 
hematomal tissue; CG, contralateral gray matter; CW, contralateral white 
matter. 



Results 

Screening of DEGs. Following data pre-processing, the gene 
expression profile data with higher normalization (Fig. 1) were 
selected for DEG analysis. The DEGs were analyzed using the 
R language limma package and BH multiple test correction. 
As a result, 399 DEGs were obtained between PH and CG, and 
756 DEGs between PH and CW according to the threshold of 
FDR<0.05 and llogFCI>l, indicating spatial specificity of the 
gene expression (Fig. 2). Furthermore, there were 35 common 
DEGs between the two groups (Fig. 3). Among the common 
DEGs, the expression change tendency for a number of the 
genes was similar [sphingosine kinase-1 (SPHKl) was down- 
regulated in PH samples compared with CG and CW]), but 
for a number was opposite [ephrin-A5 (EFNA5) and HTR2A 
were significantly upregulated in PH samples compared with 
CG but downregulated compared with C W] . 

Co-expression network. The co-expression network for the 
DEGs between the PH vs. CG and PH vs. CW groups were 
constructed utilizing String (Fig. 4). The results screened 
several interactive pairs, including inositol 1,4,5-trispho- 
sphate receptor, type 1 (ITPRl)-ryanodine receptor 2 
(RYR2) and Ca^Vcalmodulin-dependent protein kinase II a 
(CAMK2A)-MAPK1 in the PH vs. CG group, and integrin a5 
(ITGA5)-laminin Bl (LAMBl) in the PH vs. CW group. 

Enriched Gene Ontology ( GO) terms of genes in the co-expres- 
sion network. The DEGs in the co-expression network were 
enriched into GO terms using FuncAssociate. The results 
indicated that the GO terms enriched by the two groups were 
different (Table II). Compared with CG, DEGs for PH were 
significantly enriched in calcium ion transport (GO;0006816), 
ion transport (GO:0006811) and synaptic transmission 
(GO:0007268). Whereas compared with CW, the DEGs of PH 
were significantly associated with blood vessel development 




Figure 2. DEGs of the two groups. The two circles represent the DEG set of 
PH vs. CW and PH vs. CG. The red numbers indicate number of DEGs. PH, 
perihematomal tissue; CG, contralateral gray matter; CW, contralateral white 
matter; DEGs, differentially expressed genes. 



Value 



PH VS.CW 
i_ 



PH VS. CG 
1 f ' 1 




Figure 3. Hierarchical clustering of gene expression. Red represents over- 
expression and green represents low expression. The left shows PH vs. CW 
and the right PH vs. CG. PH, perihematomal tissue; CG, contralateral gray 
matter; CW, contralateral white matter. 



(GO:0001568), vasculature development (GO:0001944) and 
regulation of cell proliferation (GO:0042127). 

Pathway enriched by DEGs in network. Pathway analysis for 
DEGs of the two groups was also performed using EASE. As 
demonstrated in Table III, the four pathways were identified 
to be significantly correlated with PH when compared with 
CG, including the Ca^* signaling pathway (hsa04020), which 
was consistent with the GO enrichment analysis (Table II). A 
number of genes were shared by enriched GO terms and path- 
ways, including SLC8A2, calcium channel, voltage-dependent, 
T-type, all subunit (CACNAII), ITPRl, PRKCB, RYR2, 
CACNAIE, PPP3CA and CAMK2A. Six pathways were 
enriched for DEGs between PH and CW, among which extra- 
cellular matrix (ECM)-receptor interaction (hsa04512) was the 
most significant pathway. ECM-receptor interaction is most 
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Table II. Enriched Gene Ontology terms of differentially expressed genes in two interaction networks. 



A,PHvs.CG 


Term 


Function 


FDR 


Oenes 


OW.UUUOo 1 o 


Ctilciuni ion triinsport 




'iir'SA? r^Af^NTAiT r'Ar'NTR? r^Af^Nrn'^ r^Ar'NrA9r»'^ 

OLV^OrtZ, i 1, V^/W^iNDZ, V^rtV^iNrtZLy J , 

ITPRl, PRKCB, TPTl , RYR2, CACNAIE, PPP3CA, TRPV6, 
CAMK2A,JPH1 


OW.UUUOo 1 i 


Ion trtinsport 




KCNS 1 , SLC 1 A6, ATP5S , TPTl , TRPV6, PPP3CA, 
SLC30A3, JPHl, CAMK2A, OABRA2, SLC8A2, KCNBl, 
CACNAII, ORIA3, CACN03, CACNA2D3, ANKH, 
KCN VI, ITPRl, PRKCB, SLC17A7, KCNJ4, ORIA2, 
ORIAl , RYR2, CACNAIE 


GO:0007268 


Synaptic transmission 


5.94E-04 


DLOAPl, OABRA2, STXl A, DLOAP2, NRXN3, CACNB2, 
PAH, VIPRl, SLC17A7, MAPKl, KCNQ5, ORIA2, ORIAl, 
CHRMl , SLC1A6, CACNAIE, PPP3CA, UNC13C, HTR2A 


OW.UUU / ZD / 


Cell-cell signaling 


U .UU 1 U i / 


SLC1A6, PPP3CA, FOFl, OABRA2, DLOAPl , STXIA, 
EFNB3, NRXN3, DLOAP2, EFNB2, SLC17A7, MAPKl , 
ORIA2, ADM, ORIAl, CHRMl, EFNA5, CACNAIE, 
UNC13C,HTR2A 


GO:0019226 


Transmission of 
nerve impulse 


0.001418 


DLOAPl, OABRA2, STXIA, DLOAP2, NRXN3, CACNB2, 
PAH, VIPRl , KCNMB2, SLC17A7, MAPKl , KCNQ5, 
ORIA2, ORIAl , CHRMl , SLC1A6, CACNAIE, PPP3CA, 
UNC13C,HTR2A 


KjKJ .\J\J d\j\j\j 1 


Metal ion transport 


n nn 1 one 


C ACNA2D3 , ITPR 1 , KCNV 1 , KCNMB2 , PRKCB , 
SLC17A7, KCNQ5, KCNJ4, KCNS 1, TPTl , RYR2, 
CACNAIE, TRPV6, PPP3CA, SLC30A3, JPHl, CAMK2A 


GO:0006812 


Cation transport 


0.002788 


ATP6V0E1 , SCN3B, CACNB2, KCNMB2, KCNQ5, KCNS 1 , 
ATP5S, TPTl, SLC30A3, PPP3CA, TRPV6, JPHl, CAMK2A, 
SLC8A2, KCNBl , CACNAII, CACN03, CACNA2D3, 
KCNV 1, ITPRl, PRKCB, SLC17A7, KCNJ4, RYR2, 
CACNAIE 


00:0015674 


Di-, tri-valent 
inorganic cation 
transport 


0.002936 


SLC8A2, CACNAII, CACNB2, CACN03, CACNA2D3, 
ITPRl , PRKCB, TPTl , RYR2, CACNAIE, PPP3CA, TRPV6, 
CAMK2A, JPHl 


00:0030182 


Neuron 
differentiation 


0.036294 


DCC, NRPl, EFNB3,NRXN3, NTN4, DPYSL5, PTPRR, 
SLIT2, OPR98, ARX, NTRK3, FEZF2, DLXl , ADM, BAOl , 
ORB 1 , S 1 OOB , NEUROD2 , EFN A5 , OLFM3 


00:0022836 


Oated channel activity 


0.003996 


OABRA2, SCN3B, KCNBl, CACNAII, CACNB2, 

C ACN03 , ORIA3 , C ACNA2D3 , ITPR 1 , KCNV 1 , KCNMB2 , 

KCNQ5 , KCNJ4 , KCNS 1 , ORIA2 , ORIA 1 , RYR2 , C ACN A 1 E 


00:0005216 


Ion channel activity 


0.005017 


OABRA2, SCN3B, KCNBl, CACNAII, CACNB2, 
CACN03, ORIA3, CACNA2D3, ITPRl , KCNVl , FXYD7, 
KCNMB2, KCNQ5, KCNJ4, KCNSl, ORIA2, ORIAl, 
RYR2, CACNAIE, TRPV6 


00:0022838 


Substrate specific 
channel activity 


0.007815 


OABRA2, SCN3B, KCNBl, CACNAII, CACNB2, 
CACN03, ORIA3, CACNA2D3, ITPRl , KCNVl , FXYD7, 
KCNMB2, KCNQ5, KCNJ4, KCNSl, ORIA2, ORIAl , 
RYR2, CACNAIE, TRPV6 


00:0015267 


Channel activity 


0.01281 


OABRA2, SCN3B, KCNBl , CACNAII, CACNB2, 
CACN03, ORIA3, CACNA2D3, ITPRl , KCNVl , FXYD7, 
KCNMB2, KCNQ5, KCNJ4, KCNS 1 , ORIA2, ORIAl , 
RYR2, CACNAIE, TRPV6 


00:0022803 


Passive transmembrane 
transporter activity 


0.013258 


OABRA2, SCN3B, KCNBl , CACNAII, CACNB2, 
CACN03, ORIA3, CACNA2D3, ITPRl , KCNVl , FXYD7, 



KCNMB2, KCNQ5, KCNJ4, KCNSl, ORIA2, 
ORIAl , RYR2, CACNAIE, TRPV6 
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Table II. Continued. 



B,PH vs.CW 



Term Function FDR Genes 



GO:0001568 Blood vessel development 5.290E-04 

00:0001944 Vasculature development 6.480E-04 

00:0042127 Regulation of cell 0.001 
proliferation 

00:0006952 Defense response 0.002 

00:0008283 Cell proliferation 0.003 

00:0016310 Phosphorylation 0.003 



00:0006793 Phosphorus metabolic 0.003 

process 



00:0006796 Phosphate metabolic 0.003 

process 



00:0008285 Negative regulation of 0.004 

cell proliferation 

00:0006468 Protein amino acid 0.004 

phosphorylation 

00:0043085 Positive regulation of 0.005 

catalytic activity 

00:0032268 Regulation of cellular 0.005 

protein metabolic process 

00:0044093 Positive regulation of 0.006 

molecular function 

00:0051174 Regulation of phosphorus 0.007 

metabolic process 

00:0019220 Regulation of phosphate 0.007 

metabolic process 

00:0042981 Regulation of apoptosis 0.007 



00:0043067 Regulation of 0.007 

programmed cell death 

00:0010941 Regulation of cell death 0.008 



EDNRA, OPXl , EMCN, CTOF, HMOXl, DHCR7, COLl A2, QKI, 
ELK3 , CDH5 , SC02 , ANXA2 

EDNRA, OPXl , EMCN, CTOF, HMOXl, DHCR7, COLl A2,QKI, 
ELK3 , CDH5 , SC02 , ANXA2 

FOFR2 , DLC 1 , TESC , KAT2B , HOLS 1 , PTOS 1 , SPHK 1 , 
CDH5, ARNT, TRIB 1 , MYCN, EDNRA, AZOPl , OPXl , VDR, 
PTOES , HMOX 1 , DHCR7 , FABP3 , OPNMB , LAMB 1 , KLF4, SC02 

TF, A2M, YWHAZ, ILIRI , LYZ, HLA-C, HLA-B, PXK, STAT3, 
OCH 1 , VDAC 1 , CD 1 63 , CD83 , INHB A, HMOX 1 , IL 1 RAPE 1 , 
RAB27A, FNl , HLA-DRA, SC02 

BCATl , NANOO, MAP2K1 , NASP, ZEB2, FKBPl A, HPRTl , 
EDNRA, OPX 1 , UHRFl , DAB2 , NDE 1 , HMOX 1 , OAB 1 , CKS2 

FOFR2, IRAKI , TRPM6, MAP2K1, STK24, NEKl , P1K3C2B, 
PRKCH, STK17A, PXK, PRKX, TRIB 1 , RPS6KA3 , MAP4K5, 
PLK2, SCYL2, OAB 1 , RELN, MERTK, POKl , LIPE, UOP2 

FOFR2, NEKl , STK17A, PXK, PRKX, TRIBl , DUSP14, OABl, 
IRAKI, TRPM6, MAP2K1 , STK24, PIK3C2B, PRKCH, DUSP5, 
RPS6KA3, MAP4K5, PLK2, SCYL2, RELN, POKl , MERTK, 
LIPE,UOP2,DUSP7 

FOFR2, NEKl , STK17A, PXK, PRKX, TRIB 1 , DUSP14, OAB 1 , 
IRAKI, TRPM6, MAP2K1 , STK24, PIK3C2B, PRKCH, DUSP5, 
RPS6KA3, MAP4K5, PLK2, SCYL2, RELN, POKl , MERTK, 
LIPE,UOP2,DUSP7 

DLC 1 , TESC , K AT2B , CDH5 , TRIB 1 , VDR , AZOF 1 , PTOES , 
HMOXl, FABP3, OPNMB, KLF4, SC02 

FOFR2, IRAKI , TRPM6, MAP2K1 , STK24, NEKl , PRKCH, 
STK17A, PXK, PRKX, TRIB 1, RPS6KA3, MAP4K5, PLK2, 
SCYL2, OAB 1 , RELN, MERTK, LIPE 

DLC 1 , IRAK 1 , M AP2K 1 , SPHK 1 , ZEB2 , PMAIP 1 , TPM 1 , EDNRA , 
VDR, MAP4K5, PSMD12, OSPTl, IFT57, OABl, OCH 1, RELN 

HSP90AB 1 , DLCl , A2M, HOLS 1 , EIF5, ZEB2, FKBPl A, EDNRA, 
MAP4K5, PSMD12, TIAl , OAB 1 , QKI, HSPB 1 , PPP1R15A 

DLCl , IRAKI , MAP2K1 , SPHKl , FKBPl A, ZEB2, PMAIPl , 
TPMl , OCHl , EDNRA, VDR, MAP4K5, PSMD12, OSPTl , 
IFT57, OABl, RELN 

DLC 1 , IRAK 1 , MAP2K 1 , HOLS 1 , SPHK 1 , ZEB2 , FKBP 1 A, 
APLP2, TRIB 1 , EDNRA, INHB A, MAP4K5 , OAB 1 , CKS2, RELN 

DLCl , IRAKI , MAP2K1 , HOLS 1 , SPHKl , ZEB2, FKBPl A, 
APLP2, TRIB 1 , EDNRA, INHB A, MAP4K5 , OAB 1 , CKS2, RELN 

DLCl , IRAKI , YWHAZ, PREXl , SPHKl , ACTNl , STK17A, 
PMAIPl , OCHl , INHB A, OPXl , VDR, OSPTl , HMOXl , TIAl , 
IFT57, HSPBl, CTSB, API5, RAB27A, SC02 

DLCl , IRAKI , YWHAZ, PREXl , SPHKl , ACTNl , STK17A, 
PMAIPl , OCHl , INHB A, OPXl , VDR, OSPTl , HMOXl , TIAl , 
IFT57, HSPBl, CTSB, API5, RAB27A,SC02, DLCl 

IRAKI , YWHAZ, PREXl , SPHKl , ACTNl , STK17A, PMAIPl , 
OCHl , INHBA, OPXl , VDR, OSPTl , HMOXl , TIAl , IFT57, 
HSPB 1 , CTSB , API5 , RAB27A, SC02 



FDR, false discovery rate adjusted by BH multiple test correction method; PH , perihematomal areas; CG, contralateral grey matter; CW, contralateral 
white matter; BH, Benjamini Hochberg; GO, Gene Ontology; DEGs, differentially expressed genes; CAMK2A; Ca-Vcalmodulin-dependent protein 
kinase II a; RYR2, ryanodine receptor 2; ITPRl, inositol 1,4,5-trisphosphate receptor, type 1; LAMBl, laminin Bl; CACNAl calcium channel, 
voltage-dependent, T-type, all subunit; COLl A, collagen type 1 , a2; FNl , fibronectin 1; SPHKl , sphingosine kinase- 1; EFNA5, ephrin-A5. 
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Table III. Enriched pathways of differentially expressed genes in two interaction networks. 



A.PHvs.CG 



Pathway 


FDR 


Genes 


hsa04020: Calcium signaling pathway 


7.320E-05 


SLC8A2, ADORA2B,CACNAlI, SPHKl, 

HTR4, ITPRl , PRKCB, CHRMl , RYR2, CACNAIE, PPP3CA, 

CAMK2A, HTR2A 


hsa04010: MAPK signaling pathway 


3.020E-04 


MEF2C, NLK, CACNAII, PTPRR, CACNB2, CACNG3, 
CACNA2D3, DDIT3, PRKCB, MAPKl , 
DUSP16, CACNAIE, PPP3CA, FGFl, RAPGEF2 


hsa04720: Long-term potentiation 


0.001 


MAPKl , GRIA2, GRIAl , PPP3CA, CAMK2A, ITPRl , PRKCB 


hsa04080: Neuroactive ligand-receptor interaction 


0.007 


GABRA2, ADORA2B, GRIA2, PTH2R, THRB, GRIAl, 
CNR 1, CHRMl, HTR4, GRIA3, VIPRl , HTR2A, SLC8A2, 
CACNAII, CACNB2, CACNG3, CACNA2D3, ITPRl, PRKCB, 
TPTl , RYR2, CACNAIE, PPP3CA, TRPV6, CAMK2A, JPHl 


B.PHvs.CW 


Pathway 


FDR 


Genes 


hsa04512: ECM-receptor interaction 


0.019 


ITG A5 , COL 1 A2 , RELN , LAMB 1 , SDC2 , FN 1 


hsa05330: Allograft rejection 


0.028 


HLA-DRB 1 , HL A- A , HL A-C , HL A-B , HL A-DRA 


hsa04510: Focal adhesion 


0.031516 


M AP2K 1 , ITG A5 , COL 1 A2 , ACTN 1 , RELN , 
LAMB 1 , VASP VCL, FNl 


hsa04514: Cell adhesion molecules 


0.032 


HLA-DRB 1, HL A- A, CNTN2, HLA-C, HLA-B, 
CDH5, SDC2, HLA-DRA 


hsa05332: Graft-versus-host disease 


0.05 


HLA-DRB 1 , HL A- A , HLA-C , HLA-B , HLA-DRA 


hsa04940: Type I diabetes mellitus 


0.042 


HLA-DRB 1 , HL A- A , HLA-C , HLA-B , HLA-DRA 



FDR, false discovery rate adjusted by Benjamini-Hochberg multiple test correction method; PH, perihematomal areas; CG, contralateral grey matter; CW, 
contralateral white matter; CAMK2A; Ca'Vcalmodulin-dependent protein kinase II a; RYR2 , ryanodine receptor 2; ITPRl , inositol 1 ,4,5-trisphosphate 
receptor, type 1; LAMBl , laminin Bl; CACNAl calcium channel, voltage-dependent, T-type, all subunit; COLIA, collagen type 1, a2; FNl , fibro- 
nectin 1 ; SPHK 1 , sphingosine kinase- 1 . 



A B 




Figure 4. Interaction network of DEGs. (A) The network for PH vs. CG and (B) the network for PH vs. CW. PH, perihematomal tissue; CG, contralateral gray 
matter; CW, contralateral white matter; DEGs, differentially expressed genes. 
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commonly associated with cell proliferation and apoptosis. 
Therefore, the pathway analysis results of the DEGs in PH vs. 
CW group were also in accordance with the GO pathway [the 
common genes were collagen type 1, a2 (COL1A2), LAMBl 
and ACTNl]. 

Discussion 

The present study identified >350 genes that were signifi- 
cantly differentially expressed following ICH in PH regions 
compared with contralateral healthy tissues (grey matter and 
white matter) in human samples according to the threshold 
of FDR<0.05 and llogFCI>l. The majority of DEGs were 
further revealed to be involved in the networks and pathways 
correlated with the calcium signaling pathway, cell prolif- 
eration, regulation of apoptosis, neuron differentiation or 
ECM-receptor interaction. 

As brain metabolism is almost entirely dependent on 
the oxidation of glucose delivered by the blood, obstructed 
blood flow during ICH may cause a rapid decline in tissue 
metabolism that culminates in neuronal cell death. Although 
multiple factors are involved, disturbances in neuronal 
ionic homeostasis may be fundamental for neuronal cell 
swelling and death (24). Dissipation of ionic gradients causes 
neuronal depolarization and the release of neurotransmitters 
from intracellular stores, including glutamate. Glutamate 
provokes further depolarization due to activation of the 
a-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid 
and N-methyl-D-aspartate glutamate receptors, as well as 
voltage-gated Ca^* channels, thus triggering Ca^"^ influx 
and release from extracellular spaces. In addition, Ca^"^ is 
mobilized within the cell from intracellular Ca^* stores 
through two types of Ca~* release channels: Ryanodine 
receptors (RyRs) and inositol 1,4,5-trisphosphate receptors 
(IP3Rs) (25,26). Such Ca^"^ accumulation and excitation 
contributes to the generation of reactive oxygen and nitrogen 
species, and the release of apoptotic factors, ultimately 
causing cell death through excitotoxicity (27). As expected, 
the present results implied that CACNAII, CAMK2A, 
RYR2 and ITPRl were significantly upregulated in PH 
areas. There is also evidence that intracellular Ca^* release 
channels may cooperate, leading to positive feedback during 
activation (28). Agonist-dependent activation of IP3Rs may 
promote activation of RyRs, amplifying and shaping the 
resulting Ca'*-signal (28). The ITPR type 2 channel activity 
is modulated by Ca^VCaMKII-mediated phosphorylation of 
serine 150 of the IP3R2 (29). As expected, the present study 
also found an interaction between RYR2 and ITPRl (Fig. 4). 

Consistently several genes that regulate neuronal cell 
proliferation were downregulated but the pro-apoptotic 
genes were upregulated in the PH area. For example, SPHKl 
is an enzyme that catalyzes the phosphorylation of sphingo- 
sine to form sphingosine-l-phosphate (SIP). SIP regulates 
diverse biological processes including proliferation and 
survival following binding to SIP receptor-1 to 5 (SlPl to 
S1P5) (30,31). Furthermore, translational investigations 
suggest a profound impact of SIP administration in the 
modulation of edema formation in disease state in which 
increased vascular permeability is a hallmark feature (32). 
Isoflurane delays the development of early brain injury 



following subarachnoid hemorrhage through upregulating 
SPHKl expression (33). 

EFNA5 is reported to be upregulated in reactive astrocytes 
following stroke, which limits axonal sprouting from cortical 
neurons and motor recovery. Blockade of EFNA5 signaling 
using a unique tissue delivery system induces the formation 
of a novel pattern of axonal projections in motor, premotor 
and prefrontal circuits, and mediates neuronal recovery (34). 
As expected, EFNA5 was also upregulated in PH regions 
compared with normal CG, but it was downregulated compared 
with CW (conveys axons), indicating the ephrin/Eph A expres- 
sion gradients in the central nervous system (35). 

PH edema is directly toxic to neurons as it disrupts the 
blood-brain barrier (36). The blood-brain barrier primarily 
consists of endothelial cells with specialized tight junctions 
lining the blood vessels, astrocytic end-feet surrounding 
the blood vessels, and pericytes embedded in the basement 
membranes between the endothelial cells and the astrocytes. 
This dynamic structure is highly regulated by interac- 
tions between its cellular and ECM components (including 
collagen, laminin, fibronectin, hyaluronan and proteoglycan) 
along with integrin receptors (37,38). In the present study, the 
ECM-receptor interaction pathway was also identified and 
different ECM fibrous proteins [COL1A2, LAMBl, syndecan 
2 and fibronectin 1] and integrins (ITGA5) were downregu- 
lated. 

In conclusion, the present study provides several potential 
targets to decrease hematoma expansion and alleviate neuronal 
cell death following spontaneous intracerebral hemorrhage. 
However, further studies are required to validate this evidence 
and determine its clinical utility. 
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